Computable Integrability. 
Chapter 5: Factorization of LPDOs. 
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1 Introduction 



Let us notice first that different definitions of integrability, as a rule, use lin- 
earization of initial equation and/or expansion on some basic functions which 
are themselves solutions of some linear differential equation. Important fact 
here is that linearization of some differential equation is its simplification 
but not solving yet. For instance, in case of linear Schrddinger equation, 
i^xx + k 2 ip = uip, we are not able to find its solutions explicitly but only to 
name them Jost functions and to exploit their useful properties (see previous 
Chapters) . 

On the other hand, well-known fact is that for LODE with constant co- 
efficients operator itself can always be factorized into first-order factors and 
thus the problem is reduced to the solving of a few first-order LODEs: 



which are solvable in quadratures. 

In case of differential operators with variable coefficients factorization is 
not always possible but for the great number of operators BK-factorization 
gives factorization conditions explicitly which we are going to demonstrate 
in the next Section. 

2 BK-factorization 

Speaking generally, BK-factorization produces following result: in case of 
LPDO with characteristic polynomial having at least one distinct root, fac- 
torization is constructed algebraically for an operator of arbitrary order n 
while in case of some multiple roots of characteristic polynomial of LPDO, 
factorization is formulated in terms of Riccati equation(s). Factorization 
of LODO is always equivalent to solving some Riccati equation(s). Below 
explicit procedure for order 2 and 3 is briefly described. 



Let us outline here BK-factorization procedure pQ for the simplest case of 
bivariate LPDO of second order. Consider an operator 



d 

dx 



ip + \ip = f(x) 



2.1 LPDO of order 2 
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with smooth coefficients and seek for factorization 

A 2 = (pid x + p 2 d y + p 3 )(p4d x + p 5 d y + p 6 ). 

Let us write down the equations on pi explicitly, keeping in mind the rule 
of left composition, i.e. that d x (ad y ) = d x (a)d y + ad xy . 

Then in all cases 



O20 = PlP4 

an = P2P4 + PiPh 

O02 = P2P5 

a w = £(p 4 ) +p 3 p A +pm 

ooi = C(p 5 ) +P3P5+P2P6 

{aoo = C(pe) +p 3 Pe 



{2SysP) 



where we use the notation C = p\d x + p 2 d y . In generic case we assume that 
(after a linear change of variables if necessary) 

a 20 7^ and p\ — 1. 

Then first three equations of 2SysP, describing the highest order terms 
are equations in the variables p 2 ,P4, P5 and to find them we have to find roots 
of a quadratic polynomial 

7M0) := a 2O (0) 2 + on(0) + o 02 = 
and this leads to a linear system for p 4 , p 5 with as parameter: 



" 1 0" 




Pa 




020 




Pa 




"1 







020 


-0 1 








_a n _ 


1 









1 




_a n _ 



Thus 



Pi = 1 

P2 = -0 

Pa = a 2 o 

lP5 = O 2 Q0 + «11 



(2Pol) 



and choice of a root generates different possible factorizations of operator 
A 2 . 
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Having computed p 2 ,P4,P5 one can plug them into two next equations of 
2SysP 

a w = C(p 4 ) + p 3 P4 + PiPe 
a 01 = C(p 5 ) +P3P5+P2P6 

and get a linear system of equations in two variables p 3 , p 6 which can easily 
be solved 

0ai O + aoi - 0£a 2 o ~~ £( a 2o0 + «n) 
2a 2O + an 

(a 2O + aii)(aio - Ca 20 ) - a 20 (a i - £(a 2O + an)) 



Pe 



2a 2O + an 

(2) 



if 7*2(0) = 2a 2O + a n 7^ 0, i.e. is a simple root. At this point all 
coefficients pi,p 2 , •••,P6 have been computed and condition of factorization 

a o = £(pe) + P3P6 (3) 

takes form 



£ f 0"1O +api -£(2a2o 0+aii ) 1 _i_ 0aiu+aoi— £(2a 2 o0+aii) ^ 
|_ 2a 2 o0+an J 2a 2 o0+an 



aoo 

^ a2u(aui-£(a2()0+ail))+(a2O0+a-ll)(ai()-£a-2o) 
2a2o0+an 



(4) 



2.2 LPDO of order 3 

Now we consider an operator 

A 3 = a ikd j x d k y = a 30 d 3 x + a 21 d 2 x d y + a 12 d x d 2 y + a 03 dy 3 (5) 

i+fc<3 

+a 20 dl + a u d x d y + a 02 d 2 + a w d x + a 0l d y + a 00 . 
with smooth coefficients and look for a factorization 
^3 = (Pid x + p 2 d y + Ps)(P4dl + p 5 d x d y + p 6 dy + p 7 d x + p 8 d y + p 9 ). 
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The conditions of factorization are described by the following system: 



a 3 o 

021 

ai2 

a 03 
a 20 

an 

«02 

aio 
a i 
Ooo 



= PiPa 

= P2PA+P\Pb 

= P2P5+P1P6 

= P2P6 

= C(p 4 ) +P3P4 + P1P7 

= C(p 5 )+ P3P5 + P2P7 + PlP8 

= £(p 6 ) +P3P6 +P2P8 

= C(p 7 ) +P3P7 + P1P9 

= C(p 8 ) +P3P8+P2P9 

= C(p 9 ) +P3P9 



{SSysP) 



with C = pid x + p 2 d y . 



Once again we may assume without loss of generality that the coefficient 
of the term of highest order in d x does not vanish, and that the linear factor 
is normalized: 

O30 7^ 0, Pi = 1. 

The first four equations of 3Sys Var describing the highest order terms are 
equations in the variables P2,Pa,P5,Pg- Solving these equations requires the 
choice of a root — p<i of a certain polynomial of third degree. Once this choice 
has been made, the remaining top order coefficients P4,P5,p& are easily found. 
The top order coefficients can now be plugged into the next four equations 
of SSysP. The first three of these four equations will now be a linear system 
of equations in the variables P3,P7,ps which is easily solved. The next equa- 
tion is now a linear equation on variable pg which means that all variables 
Pi, % = 1, 9 have been found. The last two equations of SSysP will give us 
then the conditions of factorization. 

Namely, at the first step from 



030 = Pa 
021 = P2PA + Pb 
012 = P2P5 + Pa 
«03 = P2P6 
it follows that 

Vz{-p2) ■= o 3 o(-£>2) 3 + a 2 i(-p 2 ) 2 + a 12 (-p 2 ) + a 03 = 0. 
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As for the case of second order, taking p 2 = — 0, where is a root of the 
characteristic polynomial V 3 we get a linear system in P4,Ps,P6 with as 
parameter. Then again 

p 2 = -0, C = d x - 0d y , 

which leads to 

030 = P4 
a 21 = -0Pa + p 5 
012 = -0P5 + Pa 
i.e. 



1 

-0 





1 



Thus 



Pa 




a 30 




P4 


P5 




021 










_a i2 _ 




P6 



2 



Pi 


= 1 


P2 


= -0 


< J9 4 


= a 3 o 




= a 3O + a 2 i 


k P6 


= a 3O 2 + a 2 i0 + a.12 



At the second step, from 





1 

1 



C(p A ) +P3PA+P1P7 

£(ps) + P3P5 + P2P1 + pips 

£(pe) +P3P6+P2P8 



and (2 Pot) we get 



a 3 o 

«21 
Ol2 



(3Po/) 



a 20 - £a 3 o = P3O30 + Pi 

a lx - C(a 3O + a 21 ) = p 3 {a 3O + a 21 ) - 0Pi + Ps (3 Lin*) 

a 02 - C(a 3O 2 + a 21 + a 12 ) = p 3 (a 3O 2 + a 21 + a 12 ) - 0p 8 . 

As a linear system for p 3 , p 7 , p 8 this has determinant 

3a 3O 2 + 2a 2 i0 + a 12 = P 3 (0), 
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so if is a simple root the system has unique solution 

_ 2 (a 2 O - £030) + 0(on - £(^300 + «2l)) + ^02 ~ ^(a3O0 2 + ^210 + CL 12 ) _ 

PS 3a 3O 2 + 2a 2 i0 + a 12 

^20 — £«30 a 30 
P7 = — ' P3] 

3a 3O 2 + 2a 2i + a 12 3a 3O 2 + 2a 2i + a 12 

_ 0(a 2O - Ca 3 o) + an - £(a 3 o0 + Q21) a 3Q0 + ^21 

3a 3O 2 + 2a 21 + a 12 3a 3O 2 + 2a 2i + a 12 

In order to find the last coefficient p Q we use the next equation of (3Sys- 
Var), namely: 

a 10 = C(p 7 ) + p 3 p 7 + pip 9 , (3Lin**). 

At this point all coefficients Pi,i = 1, ...9 have been computed, under the 
assumption that is a simple root. 

At the third step from 

a i = C(p 8 ) + p 3 p 8 + p 2 p 9 ^ 
a 00 = C(p g ) +p 3 p 9 

all the necessary conditions for factorization can be written out. We do 
not do so here because the formulas are tedious and do not add anything 
to understanding the main idea. If the conditions are satisfied, the explicit 
factorization formulae could be written out as for the second-order operator. 
The difference is that in this case the polynomial defined by the highest order 
terms is of degree 3 and we have not one but two conditions of factorization. 



2.3 Constant coefficients 

• Obviously, in case of constant coefficients all formulae above can be 
simplified considerably and used for classical factorization problem of 
a polynomial. For instance, a bivariate second order polynomial 

W = X 2 - Y 2 + a w X + a iY + a Q0 

can be factorized into two linear polynomials, 

X 2 - Y 2 + a w X + a m Y + a 00 = (piX + p 2 Y + p 3 )(p 4 X + p 5 Y + p 6 ), 
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iff 

,2 



a m a 10 



a oo — i ^ 

In each case coefficients can be written out explicitly, for instance if 



a 2 -a 2 



a 00 — | 
then 

X 2 -F 2 + a 10 X + a ir + a 00 = (X + Y + ^^)(X -Y + ^±^). 

• As in case of order two, constant coefficients simplify all the formulae 
and reduce the problem under consideration to the classical factoriza- 
tion of a polynomial. For instance, a bivariate third order polynomial 

W = X 2 Y + XY 2 + a 20 X 2 + a u XY + a 02 Y 2 + a 10 X + a 01 Y + a 00 

can be factorized into the product of one linear and one second order 
polynomials, 

X 2 Y + XY 2 + a 20 X 2 + a u XY + a 02 Y 2 + a w X + a 01 Y + a 00 = 

= ( Pl X + p 2 Y + p 3 ) (p A X 2 + p 5 XY + p 6 Y 2 + p 7 X + p s Y + p 9 ) 
for instance, if 

ooi = oio + (a 20 + l)(a u - a 20 - a 02 ), 

a 00 = ( a ll — a 20 ~ Gk)2)[fll0 + a 2o( a ll — a 20 — a 02)}, 

and the result of factorization then is (with notation 7 = an— a 20 — a 02 ): 
X 2 Y + XY 2 + a 20 X 2 + a n XY + a 02 Y 2 + a 10 X + a 01 Y + a 00 = 
= (X + Y + -f)(XY - a 20 X + (a 20 - a u + i)Y + a 10 + a 2o7 ). 

3 Laplace transformation 
3.1 Main notions 

The most important question now is - what to do when conditions of factor- 
ization are violated? Do we still have a way to solve an equation C(ip) = 
corresponding to the initial operator? In order to answer these questions 
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let us re-write results of BK-factorization for generic case of second order 
hyperbolic operator as 



C- dd +ad +bd +c=l (d x + b)(d y + a)-ab-a x + c 

and corresponding LPDE as (d x d y + ad x + bd y + c)tpi = and introduce 
new function ip 2 = {d y + a)ip\. Our main goal now is to construct some new 
LPDE having ip 2 as a solution and to check its factorization property. If this 
new LPDE is factorizable, then its solution is written out explicitly and due 
to the invertibility of a transformation ipi — > ip 2 the formula for solution of 
initial LPDE can also be obtained immediately. 

Let us first introduce some definitions. 
Definition 3.1 Two operators of order n 



a = E a ^ d v and £*= E h ^ d v 



y 

j+k<n j+k<n 



are called equivalent operators if there exists some function / = f(x,y) such 
that 

fCi = C 2 of. 

Definition 3.2 Expressions 

a = ah + a x — c and b = ab + b y — c 
are called Laplace invariants. 

Lemma 3.3 Two hyperbolic operators C\ and C 2 of the form (J2J) are 
equivalent iff their Laplace invariants coincide pairwise. 

► Indeed, 

fd(ip) = fd x d y + faid x + fbid y + fa 

and 

£ 2 (/V) = d X dy{f)lP + dyfd^ + d x dy{lP)f + d y ^d x f + 

a2d x (f)ip + a 2 d x (il))f + b 2 d y (f)ip + b 2 d y (ip)f + fipc 2} 
i.e. fCi = C 2 o / iff 

dyf 
f 



ai = a 2 + = a 2 + (log/%; 
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bi = b 2 + ^ = b 2 + (\ogf) x ; 

ci = c 2 + + a2 d -A + 6 2 M = C2 + {logf)xy + (log/) x (log/) y + 

+a 2 (log/% + & 2 (log/%. 

Direct substitution of these expressions into formulae for Laplace invari- 
ants gives (we use notation ip = log/): 

61 = ai&i + a M - ci = (a 2 + <p v )(&2 + + (a 2 + w)*- 

C 2 - (f xy - <Px<Py - 0-2¥x ~ h^y = «2&2 + 0-2,x ~ C 2 = CL 2 . 

Analogously one can obtain &x = 6 2 and it means that for two equivalent 
hyperbolic operators their Laplace invariants do coincide.! 

► First of all, let us notice that two operators 

d = (d x + b 1 )(d y + a 1 ) + A 1 

and 

£2 = (d x + b 2 )(d y + a 2 ) +A 2 
can be transformed into some equivalent form 

Ci = (d x + h)d y + A l 

and 

L 2 = (d X + b 2 )dy + A 2 

(perhaps) by different functions /1, f 2 such that 

= /iA = C 2 o / 2 

and as was proven above, Laplace invariants of the initial operators coin- 
cide with those of the equivalent ones. Operator of the form 

C = (d x + b)d y + A 

has following Laplace invariants: A and —b y + A, i.e. A\ = A 2 and b\ y = b 2y 
for operators C\ and C 2 with the same Laplace invariants. It yields to 

b 1 = b 2 + (p(x), t 2 = (d x + bi + ip(x))d y + A x 
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with some arbitrary smooth function ip(x). Now operator L 2 differs from C x 
only by term <£>(x) which can be "killed" by one more equivalent transforma- 
tion, namely, for some function f'3 = fs(x) 

h = U'C 2 o / 3 = (d x + h + ifi(x) + (log h) x )d y + A x 
and choice (f(x) = f' z l f completes the proof. I 

Now let us rewrite initial operator (JZj) as 

£1 : d x d y + a x d x + b x d y + c x 

and notice that 

(d x + 61)^2 = a x ip x 

leads to 

( {d x + &i)V>2 = Ol^l ^ f + & l)^2 = Ipl 

\ (d y + ai)^i = 1P2 \ (d y + a x )ij; x = ^2 

and standard formula for log-derivative: 

e^dyt'^ = d y — <f y with 90 = log a, 

gives finally a new operator C 2 with corresponding LPDE 

£2(^2)= + a x - (log a x )y)(d x + b x ) - a x ip 2 = 0. (8) 

In order to check whether these two operators C x and L 2 are different, 
let us compute Laplace invariants of the new operator L 2 . 

a 1 = a 1 b 1 + a 1 , x -c 1 _ f a 2 = a x - a x>y - (logai)^ + a x>x 
b x = a x b x + - ci |6 2 = a x 

Now one can see that operators C x and £ 2 are not equivalent and operator 
L 2 is factorizable if d 2 = or 6 2 = (see example below). 

Definition 3.4 Transformation C x =>- £ 2 , i.e. 

a x = a x b x + a XjX - c x _ f a 2 = a x - a x>y - (log a x ) xy + ai iX 
61 = a x b x + - ci \b 2 = a x 

is called Laplace transformation. 
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If first new operator is also not factorizable, the procedure can be carried 
out for as many steps as necessary in order to get some factorizable operator. 
At the step N when the first factorizable operator is found, algorithm stops 
because the division on corresponding = is not possible any more. In 
fact, it is possible to write out formulae for Laplace transformation in terms 
of Laplace invariants only. 

Theorem 3.5 Let u n is one of Laplace invariants a n , b n obtained at the 
step n. Then 

U n+ i = 2u n + ()ogu n ) xy - u n -i. (9) 

► Indeed, due to Lemma 3.3 it is enough to regards sequence of operators 
of the form 

C n : d x d y + a n d x + b n d y + c n with a n = 

because it is easy to find some function f an (for instance, /< 
such that 

A n , = fan£n ° fa„ = 9 x d y + b n d y + C n . (10) 

From now on 

£n ' d x dy ~\~ b n 0y ~\~ c n 

and tilde-s are omitted for simplicity of notations. Now formulae for Laplace 
transformation take form 

1pn,y = -Cn^Pn+l, 4>n = [9 X + b n + (log C n ) x ]i) n+1 (11) 

and Eq.(jHJ) can be rewritten for the function ip n+ i as 

[d x + b n + (log Cn) x ) dyip n+ i + [c n + b n , y + (log c n ) xy )i/; n+1 = 0, 

i.e. 

c rt+l — c n — b n .y + (\ogC n ) xy , b n+ i — b n = (logCn)^, 
C-n+2 — Cn+l = + (Jog c n+l) xy , — b n+ \ = (logCn+i)^ 

and finally 

c n+2 = 2c n+ i + (log c n+1 ) xy - c n . (12) 
Notice that in this case, first Laplace invariant is 
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and obviously satisfies to Eq. (jl2j) . i.e. for the first invariant the statement of 
the theorem is proven. In order to prove it for the second invariant b n one 
has to choose another sequence of operators with b n = generated by some 
function fb n such that 



B 



n,0 



f K C n o f bn = d x d y + a n d x + 



(13) 



Notice that in order to obtain the recurrent formula for Laplace invariants, 
we used separation of variables x and y given by (fTTj) . Moreover, introduction 
of a new discrete variable n allows us to regards these equations as difference- 
differential ones. In order to deal with this sort of equations one needs a 
couple of definitions. 



Definition 3.6 An operator T acting on the infinite sequences of functions 

(..., ^_2,^-l,^0, ^1,^2, -,^n, •••) 

as 



Tijj n = ip n+1 , T Vn = tp- 



n-l 



is called shift operator. For convenience of notation sometimes infinite 
vector-function ip^ is introduced 

Ipoo = (..., -0-2, IpO, fa, -,1pn, •••) 

and matrix of operator T then has the following form: 



/ 



V 



1 

1 



0.. 

0.. 






1 















1 





\ 



(14) 



/ 



i.e. it is infinite matrix with all zero elements but the elements over main 
diagonal - they are equal to 1. 

Definition 3.7 Commutator C = [A, B] of two operators A and B is 
defined as 

C = AB - BA. 



Obviously, following properties hold true: 
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• [d x , d y ] = (cross-derivative rule), 

• [d x ,T] = [d y ,T]=0, 

• [T, a] = T o a — a o T = (T(a) — a) o T (Leibnitz rule analog). 

Let us now regard two operators corresponding Laplace transformations 
from Theorem 3.5 rewriting slightly formulae (jllj) in terms of shift operator: 

l/>n,y = -C n T4> n , 4> UtX + b n ^ n = T~ X ^ n . (15) 

The use of shift operator makes it possible to present their commutator 
C=[d y + c n T, d x + b n -T- l ] = [d y + cT, d x + b-T- 1 } 
omitting low index n. 

Lemma 3.8 Commutator 

C=[d y + cT, d. + b-T- 1 } 

is equal to zero iff 

c x = c(T(b)-b), b y = c-T- 1 (c). (16) 
► Indeed, by definition 

C = b y -c x T + cTb-bcT + T- l cT -c = -c + b y + T- l (c)-(c x -cT(b) + bc)T. 

Now, if C = 0, then -c + b y + T -1 (c) = and c x - cT(b) +bc = 0, i.e. (flfijl 
holds. 

If (fTB^l holds, then coefficients of C are equal to zero, i.e. C — 0. I 

At the end of this Section let us notice that in the original Eq.((Zj) two 
variables x and y played symmetrical role which can also be observed in 
commutation relation of Lemma 3.8 after appropriate gauge transformation: 

e q "(d y + c n T)e~ q " = d y - a n + T, e q "(d x + b n + T^e" 9 " = d x + c^T' 1 

with 

Cfi G q + , b n Qnxy Q>n Qny ■ 

It can be shown that q n satisfy the following equation 

q n , xy = e q ^- q - -e"" 1 . (17) 

The equation (fTTj) is usually called two-dimensional Toda lattice and 

plays fundamental role in the theory of Laplace transformations. 
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3.2 Truncation condition 



Definition 3.9 Truncation condition for Eg. ([9]) . namely 

(\ogu n ) xy = u n+ i - 2u n + «„_i, n = l, N + 1 
is defined by Dirichlet boundary conditions, i.e. uo = un+i = 0. 

1_ gives us Liouville equation 1 

(logMi)a^ + 2U\ = 



For example, case N 



while case | AT = 2 | yields to the system 

(logu l ) X y^u 2 -2ui + u 
(\ogu 2 ) xy = u 3 - 2u 2 + Ui 



(\ogUi) xy 

(lOg U 2 ) xy 



—2u\ + u 2 
iii — 2u 2 



These both cases are known to be integrable in quadratures. 

Case of arbitrary given iV corresponds to the system of equations with 
following matrix 



.4 



/ -2 
1 




A' 



V 



i o 

-2 1 
1 -2 








1 

1 




.. 

.. 

.. 

-2 1 



\ 



/ 



on the right hand. This matrix is some Cartan matrix and another 
choice of boundary conditions leads to Cartan matrices of other form. It 
is interesting to notice that all Cartan matrices can be constructed in this 
way. Moreover, it is proven that system of equations corresponding to each 
Cartan matrix is integrable in quadratures Most well-known source of 

Cartan matrices is semi-simple classification of algebras Lie - there is exists 
one-to-one correspondence between these matrices and semi-simple algebras 
Lie (0). 

Explicit expression for any Laplace invariant u n after n Laplace transfor- 
mations is given by following lemma which for simplicity is formulated 
for the special case Toda chain. 

iSee Ex.1 
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Lemma 3.10 Let us regard an infinite sequence of functions 



R}, n = 0,1,2, ... 



such that 



d x d y log d. 




d n = det(<9*<9>), i, j = 1, ...,n. 



Corollary 3.11 Sequence of functions 




is solution of Eq. while sequence of functions 



d 



n 



e 



d 



'n—l 



is solution of Eq. (fTTjl . 

Theorem 3.5 describes an infinite chain of equations corresponding to 
Laplace transformations and to start with this chain, we need nothing more 
then two invariants. On the other hand, many applications of this theorem 
are connected with some special problems in which different sort of finite 
chains are considered. In the next Section we will discuss two most usable 
ways to construct some finite chain of invariants and we close this Section 
with an example of 2-steps chain. 

Example 3.12 Let us regard operator 

£i = d x d y + xd x + 2, 
then its Laplace transformation gives 
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and operator C 2 is factorizable: 

£ 2 = d x d y + xd x + 1 = d x (d y + x). 
It is a simple task to write out explicitly solution ip 2 of LPDE 

C 2 {ip2) := d x {d y + x){ip 2 ) = (19) 
and afterwards solution ipi of 

:= (d x d v + xd x + 2)(i{> 1 )=0 
can be computed by formula 

fa = -d x ip2- (20) 

Indeed, introducing in Eq. ffTTJJ) notation <p — (d y + x)ip 2 we find that 
if = Y{y) is arbitrary function of one variable (...) 

i) 2 = X{x)e~ xy + I e x(y '- y) Y(y')dy'. 



3.3 Periodic closure 

Definition 3.13 Classical periodic closure for equation on Laplace 
invariants 

(log u n ) xy = u n+ i - 2u n + u n -i, n = l, N 
is denned by periodic boundary conditions, i.e. u n+ N = u n . 

In this case Cartan matrix is replaced by matrix and for > 3 
its form is 



A 



( -2 1 
1-210 
1-21 



N 



1 \ 






V 





1 



1 




-2 1 
1 -2 



(21) 



Notice 2 that matrix A^ is degenerated. It will be shown below that 
for N — 1 initial equation can be solved explicitly, while for iV = 2, 3 initial 



2 Ex.2 
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system of equations allows some reduction to one scalar equation. 



Let us regard first case | N = 1 1 , it yields to (log U\) xy = u 2 — 2u\ + m = 



and obviously u\ = gi(x)g2(y) with arbitrary smooth functions gi(x), g 2 (y). 



Case | TV = 2 | is more interesting due to the huge amount of applications 
(surfaces with constant curvature, relativity theory, etc. ) and gives rise to 
the system of equations 

0ogui) aw = 2(ui-U2), 
(log u 2 ) xy = 2{u 2 -Ulj 

with many important properties: it has conservation laws, symmetries, 
soliton-type particular solutions, etc. In particular, the reduction of this 
system can easily be constructed to the one scalar equation: 

(logui)aq, + (\ogu 2 )xy = =>- (log u x u 2 ) xy = u x u 2 = X{x)Y(y) 

with two smooth arbitrary functions X(x),Y(y). Suppose now u\u 2 = 1, 
then setting u\ = e 6 , u 2 = e~ e we find solutions of Svs. (j22j) from solutions 
of equation 

^ y + sinh^ = 0. (23) 
This equation is called sinh-Gordon equation 



Case = 3 corresponds to the system 



(logWi)a^ = U 2 + U 3 - 2ui 

(logu 2 ) xy = Ul + u 3 - 2u 2 
(logu 3 ) xy = u 1 + u 2 - 2u 3 



which allows reduction u\ = u% 

(\ogUi) xy = u 2 + ui- 2ui 
(\ogu 2 ) xy = Ui + Mi - 2u 2 
(log ui) xy = ui + u 2 - 2m 



(logWi)^ = U 2 -Ui 

(log u 2 ) xy = 2(tti - u 2 ) 



and it can be treated analogously with the case above: 
U\U 2 Ui = 1, u 2 = e , u\ — e°, u 3 — . 



6 xy = e- 2d - e y . (24) 
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Eq.([24j) is called Tzitzeica equation and its solutions give solution of 
initial system of equations. Tzitzeica equation is also very important for 
various applications [Hj. 

For both Eq. (j2Hj) and Eq. (j24j) their general solutions are not available 
as well as for the case of general N. Attempts to solve appearing systems of 
equations directly demand some tedious technique of inverse scattering and 
produce partial solutions with singularities. Method to reduce the initial sys- 
tem to one equation simplifies drastically construction of smooth solutions. 

At the end of this section let us notice that the same three equations 
which were obtained while studying truncated and periodical cases - namely, 
Liouville, sinh-Gordon and Tzizeica equations - do appear together in some 
other context 

Theorem. Nonlinear PDE of the form u xy = f(u) has higher symmetries 
iff one of three cases take place: 

f = e u , f = e u + e- u , f = e u + e~ 2u . (25) 



As we know already (see Chapter 1) that integrability of a differential 
equation is intrinsically related to its symmetry properties. Of course, pos- 
session of a symmetry does not mean that equation is integrable in some 
sense but this fact gives us a good hint on what equations might be inte- 
grable. Moreover, for some classes of differential equations it is proven [Z| 
that integrability is equivalent to some well-defined symmetry properties. 

From this point of view the theorem above justifies hypothesis that two of 
these three equations (integrability in quadratures of Liouville equation we 
have demonstrated already) have good integrability properties. As we will 
show in the next section, integrability of these PDEs can be reduced to the 
integrability of well-know ODEs. 

Thus, in contrast to the truncation condition which leads to integrability 
in quadratures for arbitrary order N, periodical closure generates a more 
complicated situation. Here we are not able to get the answer directly in 
terms of invariants but we need first to study the properties of solutions ip n 
at the all N steps simultaneously. 
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3.4 Separation of variables 



Obviously, in case of periodical closure with period N functions ipi and ipN+\ 
satisfy the same equation which was presented in Theorem 3.6: 

^n,xy + + Cn^n = 0, V 71 = 0, ±1, ±2, ... (26) 

or, in matrix form, 

(0 X dy + {b}dy + {C})^ = 

with diagonal matrices {b} and {c}. This fact allows us to regard finite 
vector-function if), 

1p = Oi, . . . IpN ) j 

which defines completely all the properties of the initial infinite system of 
equations. Notice that the fact of coincidence for coefficients of two LPDEs 
does not mean that their solutions Vv+tv and ip n also do coincide - they 
might differ, for instance, by a constant multiplier. Therefore we need now 
some notion of periodic closure for solutions. 



Definition 3.14 Bloch periodic closure is defined for the components 
of function ip as follows: 

ip n+N = k N ip n . 

Notice once more that truncation and classical periodical closure are de- 
fined for Laplace invariants, i.e. for coefficients of Eqs.()26|). while Bloch 
closure deals with solutions of the same equations. 



Corollary 3.15 

form 



Shift matrix in case of Bloch periodical closure has 



T, 



N 



/o 








k N 



where k e C is a free parameter. 



1 
1 







1 











o\ 




1 





(27) 



20 



Corollary 3.16 Basic Lemma 3.8 holds true also for the case of periodic 
closure, i.e. commutator of Laplace transformations 

C = [d y + cT N , d x + b- T/] = 

iff 

which can be checked directly. 

Example 3.17 As it was shown above, in case of N = 2 which corre- 
sponds to the classical periodic closure of u n , chain of invariants degenerates 
into sinh-Gordon equation and this closure has the form f!22D . Corollary 3.14 
allows us to construct a connection between closure for invariants and closure 
for solutions oflfTTj) : 



















-ci 








X 


. 1 - & 2. 








fa. 


y 


-k 2 c 2 







Then 

^l* + = k~ 2 ip 2 , ipi = h^2 + fax, 

and excluding ipi or ip 2 we get an equation of the second order on one scalar 
function, for instance 

*p2,xx + Oi + b 2 )ip2x + {hx + &1&2 - k~ 2 )ip 2 = 0. 
This equation is obviously equivalent to linear Shrodinger equation 

ipxx = {\ + u)ip, A = k~ 2 , 

i.e. Sinh-Gordon equation is S-integrable. This important fact plays role, 
for instance, while constructing surfaces of constant curvatures (see very ex- 
haustive review jH]). 



4 General invariants and semi-invariants 

Before discussing the notion of general invariant, let us notice that arbitrary 
LPDO of second order, A2, can be represented in the form of factorization 
with reminder 

M = (pid x + p 2 d y + p 3 )(p4d x + p 5 d y + p 6 ) - l 2 
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where reminder l 2 is defined by (J3J): 



0a w + a 01 - £(2a 2 o0 + an) 



0a w + a i - £(2a 2O + a n ) 



x 



2a 2 o0 + an 



2a 2O + an 



x a2o(aoi ~ £(a2Q0 + a n)) + ( Q 2Q0 + Qii)(aio ~ ^20) ^ 

2a 2O + an 

Similar to the case of order two, arbitrary LPDO of third order A3 can 
be represented in the following form 

M = (pid x +p 2 d y +p 3 ) {pid 2 x + p 5 d x d y + p 6 d y + p 7 d x + p$d y +pg)- hd y - 1 31 . 

In contrast to the second order LPDO, in this case factorization with re- 
minder gives us not a function but a linear first order operator and it is con- 
venient for our further investigations to regard in this case two "reminders" 
£3 and /31 which are defined by Sys.(JEJ): 

h = a i - (pid x +p 2 d y + p 3 )p 8 - P2P9 

hi = a o - (Pid x + p 2 d y + p 3 )P9- 

In the next Section it will be shown that "reminders" l 2 , I3, Z31 are invari- 
ants of the corresponding operators uner the equivalence transformations. 

4.1 Construction of invariants 

Let us first recollect definition of two equivalent operators. 

Definition 4.1 Two operators of order n 



are called equivalent operators if there exists some function / = f(x,y) such 
that 



The definition is given for an operator £ of arbitrary order n and obviously 
any factorization 



a = E a ^ d v and ^ = E h *fi$ 




fCi = C 2 of. 



C = dC 2 



can be written out in equivalent form 



r'c o / = f- 1 Ci£ 2 o / = (r^i/Xr 1 ^ o 
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as well as sum of operators C = C\ + C 2 . 

r x c o / = /-^a + c 2 ) o f = (f- 1 ^ o /) + (/-% o /). 

Below we will take function / in a form f = e v for convenience. In order to 
formulate theorem on invariants we need following notations: 

-42a = a 2 odl + a u d x d y + a 02 d 2 + a w d x + a 01 d y + a o, 

A 2p = (pid x + p 2 d y + p 3 )(p 4 d x + p 5 d y + p 6 ) - l 2 , 

A 3a = a 30 dl + a 2 id 2 x d y + a 12 d x d 2 + a 03 <9j + a 20 d 2 x + 
+a n d x d y + a 02 d 2 + a 10 d x + a 01 d y + a 00 , 

A 3p = (Pid x + p 2 d y + p 3 ){p±d 2 x + p 5 d x d y + p 6 d 2 + 
+p 7 d x + p 8 dy + pg) - l 3 d y - hi, 

Pi = i~ x Vi ° f, Oij = f~ la ij ° f, 

Ai = j~ x A, t of, i = 2a, 2p, 3a, 3p. 

Above A 2p = A 2a , i.e. A 2p and A 2a are different forms of the same operator 
- its initial form and its form after the factorization with reminder. The same 
keeps true for A 3p and A 3a , i.e. A 3p = A 3a . 



Theorem 4.2 For an operator of order 2, its reminder l 2 is its invariant 
under the equivalence transformation, i.e. 

l 2 — l 2 - 

For an operator of order 3, its reminder l 3 is its invariant, i.e. 

^3 — ^3; 

while reminder l 31 changes its form as follows: 

hi — hi + hfy 
► Indeed, for operator of order 2 

-42a — A 2p — l 2 , 
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i.e. 

A 2a = f- 1 A 2a of = r\A 2p -l 2 )of = 

= A 2p - r\i 2 )o f = a 2p -i 2 . 

For operator of order 3 

A 3a = A 3p — l 3 d y — Z31, 

i.e. 

As a = f^Ma O f = f-\A 3p - hdy - hi) O / = 

= A 3p - f~\hi + hdy) of = A 3p - hdy - hy y - h\M 

Corollary 4.3: If h = 0, then hi becomes invariant. 

That is the reason why we call l 31 further semi-invariant. 

Corollary 4.4: If l 3 7^ 0, it is always possible to choose some function 
/ : [31 = l 3 <p y + hi = 0. 

Notice that for second order operator, if its invariant l 2 = then opera- 
tor is factorizable while for third order operator two its invariants have to be 
equal to zero, l 3 = hi = 0. On the other hand, if operator of third order is not 
factorizable we can always regard it as an operator with only one non-zero 
invariant. Of course, all this is true for each distinct root of characteristic 
polynomial, so that one expression, say, for l 3 will generate three invari- 
ants in case of three distinct roots of corresponding polynomial. Expressions 
for invariants l 2 and l 3 and also for semi-invariant l 3 \ can be easily written 
out explicitly using formulae given by BK-factorization (Section 2.1 and 2.2). 

As it was show already, for an important particular case - hyperbolic 
operator of second order in the form 

d x d y + ad x + bd y + c (29) 

- there exist two Laplace invariants which coincide pairwise for equivalent 
operators (Lemma 3.3). After rewriting hyperbolic operator in the form 

d 2 x - d 2 y + ad x + bd y + c (30) 

by appropriate change of variables, we can construct Laplace invariants as 
a simple particular case from the formulae for general invariants (see next 
Section). 
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Corollary 4.5: Two hyperbolic second order operators having the same 
normal form, ()29|) or ([30)1 . are equivalent iff their general invariants coincide. 

4.2 Hierarchy of invariants 

As it was shown above, every general invariant is a function of a distinct 
root of the characteristic polynomial and each distinct root provides one 
invariant. It means that for operator of order n we can get no more than 
n different invariants. Recollecting that BK-factorization in this case gives 
us one first order operator and one operator of order n — 1, let us put now 
following question: are general invariants of operator of order n — 1 also in- 
variants of corresponding operator of order n? 

Let regard, for instance, operator of order 3: 
-43a = A 3p = A 1 A 2a -l 3 d y -l 31 = A 1 (A 2p -l2)-hd y -h 1 = A 1 A 2p -l 2 Ai~l 3 d y -l 31 
and obviously 

A 3a = A x A 2p - l 2 A x - l 3 d y - l 3 i = A x A 2p - l 2 A x - l 3 d y - f 31 , 

i.e. l 2 is also invariant of operator A 3a . Let us notice that general invariant 
l 3 = l 3 (0^ 3 ') is a function of a distinct root 0^ of the polynomial 

V 3 (z) = a 30 z 3 + a 21 z 2 + a 12 z + a 03 

while general invariant l 2 = l 2 (0^) is a function of a distinct root 0^ of the 
polynomial 

K 2 {z) = p^z 2 +p 5 z + p 6 

with p 4 , p 5 , p 6 given by (3Pol) for = 0^ 3 \ In case of all distinct roots of both 
polynomials V 3 {z) and TZ 2 (z), one will get maximal number of invariants, 
namely 6 general invariants. Repeating the procedure for an operator of 
order n, we get maximally n\ general invariants. In this way for operator of 
arbitrary order n we can construct the hierarchy of its general invariants 

Int In— It ■■•1 ^2 

and their explicit form is given by BK-factorization. 

For instance, let us regard a third order hyperbolic operator in the form 

C = a 30 d 3 + a 2 id 2 £ d y + a\ 2 d x d 2 + a 03 dy 3 + terms of lower order (31) 
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with constant high order coefficients, i.e. = const V i + j = 3 and all 
roots of characteristic polynomial 

a 3 o0 3 + a 2 i0 2 + ai20 + a 03 = V 3 (0) 

are distinct and real. Then we can construct three simple independent general 
invariants in following way. Notice first that in this case high terms of (|31|) 
can be written in the form 

(a x d x + f3id y )(a 2 d x + (3 2 d y )(a 3 d x + (3 3 d y ) 

for all non-proportional oij,(3i and after appropriate change of variables this 
expression can easily be reduced to 

d x d y (d x + d y ). 

Let us introduce notations 

di = d x , d 2 = d y , d 3 = di + d 2 = d t , 
then all terms of the third and second order can be written out as 

C ijk = {di + CLi){dj + aj)(d k + a k ) = didjd k + a k didj + ajdid k + aidjd k + 
+(dj + aj)(a k )di + (di + a i )(a fc )(9 i + (<% + a i )(a i )9 fc + (di + a i )(9 7 - + a j )(a k ) 
with 

a 20 = a 2i a 02 = a lj a ll = 01+02 + 03 

and Cij k = C — Cij k is an operator of the first order which can be written out 
explicitly. As it was shown above, coefficients of Cij k in front of first deriva- 
tives are invariants and therefore, any linear combination of invariants is 
an invariant itself. These invariants have the form: 

for C123 we have 

a 2 a 3 + aia 2 + d y (a 3 ) + d x (a 2 ) - a w , aia 3 + aia 2 + d x (a 3 ) + d x (a 2 ) - a 01 ; 
for C312 we have 

a 2 a 3 + a x a 2 + d t (a 2 ) + d x (a 2 ) - a 10 , O1O3 + oia 2 + d t (ai) + ^(02) - a 01 ; 
for C231 we have 

a 2 a 3 + aia 2 + d y (a 3 ) + d y (ai) — a w , aia 3 + a\a 2 + d t (a\) + d y (a\) — a i. 

Direct calculation gives us three simplest general invariants of the initial 
operator C: 

hi — a 2,x — ^32 = ^3,2/ — a 2,ti hi = a 3,x ~ CLl,t- 
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Proposition 4.6 General invariants Z21, h2, ^31 are all equal to zero iff 
operator C is equivalent to an operator 

L = 6,0283 + Mi + b 2 2 + c, (32) 
i.e. 3 function / : / -1 C a f — L. 
► Obviously 

f-\d x O y O t ) o / = (0 X + (log/).) (0 y + (log /)„)($ + (log/)*) 

for any smooth function /. Notice that it is the form of an operator Cijk and 
introduce a function / such that 

0l = (l0g/)ai, 2 = (l°g/)v> 03 = (l0g/)f 

This system of equations on / is over-determined and it has solution / 

iff 

0-2,x — Q>l,y = 0, d3 >y — d2,t = 0, a^ yX — (Xi y t = 0, 

i.e. l 2 i = I32 = hi = 0. I 

■< Indeed, if C is equivalent to (|52^l . then a 2 o = a 02 = an = and obvi- 
ously hi = l m = hi = 0. 1 

At the end of this section let us notice that the reasoning above can 
be carried out for the hyperbolic operator of order n with constant leading 
coefficients and analog of Proposition 4.6 keeps true with n order terms as 

<9i<9 2 <9 3 ....<9 n 

with 

Qi = d x , d 2 = d y , d 3 = d x + d y , d A = a 4 d x + (3 4 d y , ... d n = a n d x + (3 n d y 
and terms of order n and n — 1 can be written as 

Qi,i 2 ,...,i n = {di + a t )(d 2 + a 2 )...(d n + a n ) 
while corresponding linear general invariants will take form 

kj = di(aj) - dj(ai). 
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5 Examples 

5.1 Operator of order 2 

Let us regard first for simplicity LPDO of the second order with constant 
leading coefficients, i.e. 

dij = const V(i +j) — 2 

and all roots of characteristic polynomial are distinct. Then obviously any 
root also does not depend on x, y and expressions for pij can be simplified 
substantially. Let us introduce notations 

a OO 0° = V o (0), 
a-100 1 + a i = Vi(0), 
a 2 o0 2 + an0 + a 02 = V 2 (0) 
and notice that now and V 2 {0) ^ are constants. 

Using formulae from Section 2.1 we get 



Pi = l 

J5 2 = -0 
Pl(0) 

P3 ^(0) 

P5 = a 2 o0 + On 

(a 2 o0 + an)aio — a 2 oaoi 



P6 ^(0) 
which yields to 



rl \ Q-200 + a ll \ a 20 \ 

£(p 6 ) = -^-^axo) - p^Ki) 



and invariant / 2 takes form 



; ri \ i 0200 + On a 20 , , 

H 2 {0) H 2 (0) 

. V\{0) (a 2O + a n )a 10 - a 20 a i , . 

tW) a °°- (33) 
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• Let us regard hyperbolic operator in the form 

d xx - d yy + a 10 d x + a m d y + a 00 , (34) 

i.e. a 2 o = 1, an = 0, a 2 = —1 and = ±1, C = d x — 0d y . Then Z 2 takes 
form 

_ r( 0aiQ - a 01 0a w - a m 0a w + a i 

h ~ a °° ~ A 20 } 20" 2^~~ 

which yields, for instance for the root = 1, to 

2 2 2 2 

; _ r / a W~ Ooix a 10 — a 01 _ / o a N/ a 10 — a 01x a 10 — a 01 

h - a o-L,{ ) a oo -{0 x -0 y ){ ) 

and after obvious change of variables in (jH4j) we get finally first Laplace 
invariant a 

l 2 = c — d x a — ab = a, 

where 

Qio — clqi aio + aoi 
a = , b = , c = a 00 . 

Choice of the second root, = — 1, gives us the second Laplace invariant 
b, i.e. Laplace invariants are particular cases of the general invariant 
so that each Laplace invariant corresponds to a special choice of 0. 

• Let us proceed analogously with an elliptic operator 

d xx + d yy + a>iod x + a m d y + aoo, (35) 
then = ±i, C = d x — 0d y and 

; , -r- -a \ / ±a W + O-qiI .«io + a oi 
k = aoo + {dx =F ) + i ^ 

where choice of upper signs corresponds to the choice of the root = % 
and choice of lower signs corresponds to = —%. 

5.2 Operator of order 3 

Now let us regard LPDO of the third order with constant leading coefficients, 
i.e. 

ciij = const V(z + j) = 3 
with at least one root distinct of characteristic polynomial 

O3O0 3 + Gt210 2 + Ol20 + Ct03 = V^ifi) 
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and notice that now and V' 3 {0) ^ are constants. 
Using formulae from Section 2.2 we get 



P3 



Pi = 1 
P2 = -0 

Pi = O30 



Pi 

0a ; 

Ps 



P5 = «3O0 + a 2 i 

P6 = Ct3O0 2 + ^210 + ai2 
^20 ^30 ^2(0) 



^(0) v'MV'M 

0a 2 o + an _ 0^30 + «2i ^2(0) 
^3(0) ^3(0) ^3(0) 



^(^20) Q30 ^(^2(0)) ^2(0) , Q20 «30 ^2(0) ^ 

P9_ai ° n(0) ^(0) ^(0) ^3(0)^3(0) WW 
and / 3 , Z31 are 

-(piA + pafy + p 3 )ps ~ P2P9 + a 01 = h, 

-(pid x + p 2 d y + p 3 )p9 + a 00 = hi- 

The formulae are still complicated and in order to show the use of them let 
us regard here one simple example of an operator 

B = d 2 x d y + d x d 2 y + a u d x d y + a 1Q d x + a 01 d y + a 00 , (36) 
with a 30 = a 03 = a 2 o = a 02 = 0, a 2 i = a i2 = 1. Then its invariant 

^3 = <9x«ii — «oi 

and semi-invariant 

^31 = &x a \Q ~~ a 00 

have very simple forms and gives us immediately a lot of information about 
the properties of operators of the form 1)36)1 . for instance, these operators are 
factorizable, i.e. has zero invariants l 3 = Z 3 i = 0, iff 

an = J a m dx + fi(y), a w = J a 00 dx + f 2 {y) 

with two arbitrary functions on y, f\(y) and f 2 (y). Another interesting fact is 
that if coefficient an = a n (j/) is function of one variable y, then aoo is general 
invariant and there definitely should be some nice geometrical interpretation 
here, etc. 
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6 Summary 



At the end of this Chapter we would like to notice following very interest- 
ing fact - beginning with operator of order 4, maximal number of general 
invariants is bigger then number of coefficients of a given operator, 



It means that general invariants are dependent on each other and it will be a 
challenging task to extract the subset of independent general invariants, i.e. 
basis in the finite space of general invariants. 

As to semi-invariants, notice that an operator of arbitrary order n can 
always be rewritten in the form of factorization with reminder of the form 



and exact expressions for all U are provided by BK-factorization procedure. 
The same reasoning as above will show immediately that l n is always general 
invariant, and each l n ,k-i is io-th semi-invariant, i.e. it becomes invariant in 

case if l n ,k-i = 0, Vi < 

In this paper, explicit formulae for I2, h, hi are given. Formulae for 
higher order operators can be obtained by pure algebraic procedure described 
in pQ but they are too tedious to be derived by hand, i.e. programm package 
for symbolical computations is needed. 

Already in the case of three variables, the factorization problem of a cor- 
responding operator and also constructing of its invariants becomes more 
complicated, even for constant coefficients. The reason of it is that in bivari- 
ate case we needed just to factorize leading term polynomial which is always 
possible over C. It is not the case for more then 2 independent variables 
where a counter-example is easily to find (see Ex.5), i.e. there exist some 
non-trivial conditions to be found for factorization of polynomials in more 
then two variables. 
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l n d k x + l n ,\d k x 1 + ... + l n>k -i, k<n 



7 Exercises for Chapter 5 



1. Using this particular solution u 



(x+y) 2 



of the Liouville equation 



(\ogu) xy + 2u = 
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find general solution of this equation. 

Hint: Liouville equation is invariant under the following change of vari- 
ables: 

x = X(x), y = Y(y), «(i,y) = jfpg. 

2. Prove that matrix (J2T|) is degenerate. 

3. Let Laplace invariants are equal, i.e. a — b. Prove that initial operator 

d x d y + ad x + bd y + c 
is equivalent to the operator 

d x d y + c. 

4. Let in Lemma 3.8 function w(x,y) is chosen as 

w(x,y)=X 1 (x)Y 1 (y)+X 2 (x)Y 2 (y). 

Prove that ofa = 0. 

5. Check that 

x +y +z — 3xyz = (x + y + z) [x +y +z — xy — xz — zy) 
and prove that x 3 + y 3 + z 3 is not divisible by a linear polynomial 

ax + j3y + 7Z + 6 
for any complex coefficients a, j3, 7, 5. 
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